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In a magnetization vs. temperature (M vs. T) experiment, the blocking region of a magnetic nanoparticle 
(MNP) assembly is the interval of T values were the system begins to respond to an applied magnetic field 
H when heating the sample from the lower reachable temperature. The location of this region is determined 
by the anisotropy energy barrier depending on the applied field H, the volume V, the magnetic anisotropy 
constant K of the MNPs and the observing time of the technique. In the general case of a polysized sample, 
a representative blocking temperature value Tb can be estimated from ZFC-FC experiments as a way to 
determine the effective anisotropy constant. 

In this work, a numerical solved Stoner-Wolfharth two level model with thermal agitation is used to sim¬ 
ulate ZFC-FC curves of monosized and polysized samples and to determine the best method for obtaining a 
representative Tb value of poly sized samples. The results corroborate a technique based on the T derivative 
of the difference between ZFC and FC curves proposed by Micha et al (the good) and demonstrate its relation 
with two alternative methods: the ZFC maximum (the bad) and inflection point (the ugly). The derivative 
method is then applied to experimental data, obtaining the Tb distribution of a polysized Fe^O/^ MNP sample 
suspended in hexane with an excellent agreement with TEM characterization. 


I. INTRODUCTION 

Magnetic nanoparticles (MNPs) are been extensively 
studied due to their multiple applications in technology^ 
and biomedicine^’^. Particles with sizes in the range 
[5,100] present a magnetic behaviour determined by 
its volume, shape and composition, matrix viscosity and 
temperature, among other factors. In the simplest (how¬ 
ever very useful) model, the MNPs of volume V and satu¬ 
ration magnetization Mg are considered as almost spher¬ 
ical ellipsoids with a permanent moment m = MgV and 
a preferential magnetization axis (easy axis) in which the 
anisotropy energy Fk = KV sin^ [^] is minimum, being 
K the effective anisotropy density constant and 5 the an¬ 
gle between m and the easy axis. If the MNPs are fixed 
in the matrix and separated one from each other by a dis¬ 
tance d > dipolar interactions can be neglected^ 

and the energy of the system can be expressed as the 
sum of the anisotropy energy and the Zeeman energy 
Eh = —mH cos[6>]: 

E = Ek + Eh, (1) 

with 0 the angle between m and H (fig. 1). This config¬ 
uration is usually called St oner-Wolfharth system since 
the publication of a work^ in which the authors perform 
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a numerical calculation of the M vs. H curves of ordered 
systems with different orientations i. e. systems of identi¬ 
cal MNPs with a single value of 0, and the M ys.H curve 
of a disordered system i.e. with a uniform distribution of 
(f values. Since no thermal agitation was considered by 
Stoner and Wolfharth, their calculations were made just 
finding the positions Oi of the minima of equation 1 for 
each value of H. 



FIG. I. MNP model. The energy is determined by the angle 6 
between the magnetization M and the field H, and the angle 
0 between M and the easy axis K. For calculation simplicity 
the angle (f = 0 S between K and H is used. 

In order to calculate the temperature dependence of 
the magnetic response for MNPs systems, it is neces¬ 
sary to consider the effect of thermal fluctuations that 
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allow transitions between stable configurations. Doing 
so, it is possible to simulate M vs. T experiments as the 
extensively performed Zero Field Cooling-Field Cooling 
(ZFC-FC) routine. In this kind of experiments, a sample 
is cooled from a temperature where all particles show su- 
perparamagnetic behaviour to the lowest reachable tem¬ 
perature (usually around 3K), then, a small constant 
field usually lower than SkA/m is applied, and the sam¬ 
ple is heated to a temperature high enough to observe an 
initial growth and subsequent decrease in magnetization, 
i.e. were the sample show again superparamagnetic be¬ 
haviour .The sample is then cooled again to the lowest 
temperature with the constant field still applied. 

In the ideal case of a monosized, non interacting MNPs 
sample; a narrow temperature region should exist in 
which the system performs a transition between irre¬ 
versible and reversible regimes. When heating with 
applied field, the thermal energy kT is initially much 
smaller than the anisotropy barrier KV so the magne¬ 
tization remains null. Due to the exponential depen¬ 
dence of the Neel relaxation time with temperature^, 
when kT ^ KV^ the magnetization grows rapidly until 
its thermodynamic equilibrium value, defining the afore¬ 
mentioned transition region. The Blocking Temperature 
Tb can be considered as the inflection point of this grow¬ 
ing and its experimental determination is an important 
goal of the MNPs characterization. 

Real samples always present a size dispersion, usually 
reasonably well described by a log-normal distribution. 
Different particle size implies different anisotropy bar¬ 
rier KV and therefore a different Tb for each size frac¬ 
tion, so, in real ZFC-FC experiments, the blocking re¬ 
gion is wide and the representative Tb value is not well 
defined. There are several different criteria used to de¬ 
fine a representative Tb from ZFC-FC data of polysized 
samples. Some authors maintain the inflexion point (IP) 
criterion^ while others report the maximum ZFC mag¬ 
netization temperature (MAX)^^^, being all this criteria 
still in discussion^ In an alternative approach, Micha 
et al}‘^ propose a method in which the Tb distribution is 
obtained from the T derivative of the difference between 
ZFC and FC curves. An approximated theoretical jus¬ 
tification for this method was presented by Mamiya et 

In this work, a SW model with thermal agitation is 
applied to obtain the temporal dependence of the mag¬ 
netization M(t) for an ordered system of identical MNPs 
in a similar way to previous works of Lu^^, Usov^^ 
and CarreyTemperature dependence dM{T)/dT is 
then obtained in order to numerically simulate the ZFC- 
FC curves. In contrast to the method implemented by 
Usov^^ where a stair-step approximation for the time evo¬ 
lution of the temperature was used, we consider a con¬ 
tinuous time evolution. Finally, an ordered polysize sys¬ 
tem response is simulated by linear combination of the 
monosize curves weighted by a discrete log-normal dis¬ 
tribution. 

The validity of the method proposed by Micha et al 


was verified by comparing the T derivative of this ZFC- 
FC curve with the Tb distribution obtained from the 
inflection points of each volume of the distribution. The 
resultant mean blocking temperature value < Tg > is 
then compared, for several volume distributions, with the 
commonly used criteria for a representative T^: the in¬ 
flection point temperature IP and the maximum MAX of 
the ZFC curve. 

Additionally, Micha’s method is tested with experi¬ 
mental data of a magnetite MNPs frozen ferrofluid sus¬ 
pended in hexane comparing the obtained Tb distribu¬ 
tion with the one calculated from the TEM size informa¬ 
tion. In order to obtain an ordered system, the ferrofluid 
was frozen while a large constant was field applied. 


II. MODEL 

A SW-like model with thermal agitation and zero 
width energy minima approximation was developed in 
order to obtain ZFC-FC curves of fixed MNPs with size 
dispersion. Only the simplest case of an ordered system 
was considered, with all the MNPs oriented (easy axis 
orientation) in the direction of the field. This situation 
can be achieved experimentally by freezing a ferrofluid 
sample under a sufficiently strong applied field (^ IT). 


A. Magnetization vs. time equation 

For a system of identical, fixed, non interacting MNPs 
of volume V, anisotropy constant K and saturation mag¬ 
netization Mg, with their anisotropy axes parallel to an 
external field the energy can be expressed as the sum 
of the anisotropy energy and the Zeeman energy 

E = Ek + Eh = KVSin{ef - iioM,VCos{e) 

= KV (^Sin{ef - 2hCos{e '^, 

being h = H/Hk and = \J.oMsl2K. 

In the range ^ = [0, 27r], this energy landscape presents 
two minima, of E{0) = —2KVh and E{'k) = 2KVh and a 

maximum of E(arccos(—h)) = KV (fig. 2). 

The frequency of thermal inversions between minima i 
and j is the inverse of the Neel relaxation time^^^^,: 

fij = (3) 

with /o the “intrinsic frequency”, times the Boltzman 
“success probability” depending on the ratio between 
thermal energy kT and barrier height Dij. The barrier 
between minima is symmetric for h = 0 with D^u = 
Dud = KV (naming u and d to ^ = 0 and 0 = ir direc¬ 
tions respectively) and smaller for inversion to the field 
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FIG. 2. Energy landscapes as a function of theta for a fixed 
MNP under a magnetic field at ^ = 0. Each color stands for a 
different value of the reduced field factor 2h. Blue: 2/i = 0.1, 
purple: 2h = 0.5, gold: 2h = 2. 


direction otherwise: 


/S^ud = ^(arccos(—/i)) — E^tt) 
Adu = E{diTCCOs{—h)) — E{0) 


KV 1 



KV 1 



( 4 ) 



EIG. 3. Simulation result for an ordered assembly of MNPs. 
The system is first cooled with zero field applied from a high 
temperature where all particles show superparamagnetic be¬ 
haviour (ZEG, no showed), then, the field is turned on and the 
system is heated beyond the blocking region (ZEGW). Main¬ 
taining the applied field, the system is cooled (EG). The final 
heating can be performed with (EGW) or without applied 
field (TRM). 


B. Magnetization vs. Temperature equation: ZFC-FC 
simulation 


It is a good approximation to consider the same /o 
value for both frequencies^^. 

Sample magnetization M in the direction of the ap¬ 
plied field can be expressed in terms of saturation magne¬ 
tization Ms and the number of particles per unit volume 
magnetized in each direction Nu and Nd'. 

Na^ Ms/N = M, - 1^, (5) 

with N the total number of particles per unit volume. So 
the time derivative of the magnetization can be written 
in terms of the population variation which is equal to the 
actual population times the inversion probability to each 
direction 


dM ^ ^Ms dNu 
dt N dt 


dNu 

dt 


dNd 

dt 


fd^uEd fu^dEu- 


( 7 ) 


so the time derivative of the relative magnetization m is 
determined by the transcendental equation 


= 2/oe ) {sinh(2C h) — m cosh(2C h)}. (8) 

where C = KV/kT. 


Temperature dependence of the magnetization can be 
obtained from 8 via the equation 

dm dm dt 

~dT ^ hdM' ^ ^ 

For a linear temperature variation T{t) = 5t + To, the 
magnetization derivative is 


^ ^ ^ cosh(2C/i)} 

( 10 ) 

Solving this equation by numerical methods it is pos¬ 
sible to simulate a ZFC-FC experiment for a monosize 
sample. A Matlab script based on the ODE15s^^ func¬ 
tion was developed. An example of the result for a mono¬ 
size assembly of ordered MNPs is shown in figure 3. Line 
colours stand for different parts of the routine. 

During the warming after zero field cooling (ZFCW 
for this chapter, usually called just ZFC), the exponen¬ 
tial dependence of the inversion frequency with tempera¬ 
ture in equation 3 determines a narrow “blocking region” 
wherein the MNPs, which were “blocked” at low temper¬ 
ature begin to respond to the field. Magnetization grows 
with temperature since the applied field has decreased 
the energy barrier for ^ = tt to ^ = 0 inversion. The 
magnetization increasing reverts when thermal energy is 
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FIG. 4. Comparison between ZFC-FC simulations for small 
(a = 10“®) and large (a = 0.5) dispersion systems. 



FIG. 5. Scheme of the method verification. Monosize ZFC- 
FC curves are simulated from the dM(T)/dT equations of the 
model. In one path, Tb for every particle size is calculated 
as the IP of the monosize ZFC curve. In the other path, 
a polysize ZFC-FC curve is simulated by linear addition of 
the monosize values. Then the T derivative of the difference 
ZFC-FC is calculated. 


much higher than the barrier, so the difference between 
inversion frequencies in each direction tends to disappear. 
The blocking temperature Tg of the system is then de¬ 
fined as the inflection point of the magnetization growing 
when heating. 

When the system is cooled again (FC), magnetiza¬ 
tion grows monotonically while the barrier height differ¬ 
ence between inversions becomes increasingly significant 
against thermal energy. This growth stops when thermal 
energy becomes too low for inversions to occur within 
the experimental window time. If the system is then 
heated maintaining the applied field (FCW), magnetiza¬ 
tion values are the same than FC except for the blocking 
region where there is a small increase due to the assem¬ 
bly getting closer to the equilibrium state. If the final 
warming is done with no applied field (Thermal Rema¬ 
nent Magnetism, TRM), magnetization drops to zero in 
the blocking region when thermal energy is enough for 
the wells populations to equilibrate. 

The magnetization values Mp for a polysized sample 
are obtained by linear addition of the Myi values for 
each contemplated size Vi , weighted by the corresponding 
volume and log-normal distribution LnN(Vi) value: 


Mp{T) 


Y.t^Mvi{T)ViLnN{V) 


( 11 ) 


The Vi LnN (Vi) product stands for the relative volume 
distribution. 

Figure 4 shows the comparison between ZFC-FC sim¬ 
ulations for samples with different size dispersion ex¬ 
pressed as the scale parameter a of the log-normal num¬ 
ber distribution. A much wider transition region can be 
seen for the bigger dispersion so the different aforemen¬ 
tioned criteria would define very separated values for a 
representative Tg. 


III. BLOCKING TEMPERATURE DETERMINATION 
A. Micha’s method verification 

In order to verify Micha’s method, several polysize 
ZFC-FC experiments were simulated using differents pa¬ 
rameter sets varying a and the mean radius. For each 
one of the used sets, the T derivative of the ZFC-FC dif¬ 
ference was calculated. Then, the Tb distribution was 
obtained from the monosize curves that were added to 
construct the polyzise simulation in equation 11: a ZFC 
curve was calculated for each class of the size distribution 
so each Tb class comes from a volume class, maintain¬ 
ing the same relative height. Also IP and MAX values 
of the polysize ZFC curve were calculated and compared 
with the mean value < Tg > of the distribution in each 
simulation (fig. 5). 

In all cases, the Tb distribution and the ZFC-FC 
derivative are identical. Figure 6 shows the results for the 
simulation with 4.5 nm mean radius, a = 0.5, 16k J/m^ 
anisotropy constant and a AK/min heating rate. 

Also, for a set of ZFC-FC curves calculated with the 
same mean volume, saturation magnetization, heating 
rate and anisotropy constant, by increasing scale param¬ 
eter cr, <Tb > stays constant while the polysize curve IP 
shifts to smaller temperatures and MAX shifts in oppo¬ 
site direction. Figure 7 shows the results for 4.5nm mean 
radius, \6kJjw? anisotropy constant, ^Kjmin heating 
rate and a = [0.1, 0.6]. 

This behaviour is the same in the hole studied size 
range. By normalizing IP values by the Tb mean, all 
points fall in the same curve as shown in figure 8 while 
the variation for MAX is small. 

Varying the heating rate and K does not affect the 
relation IP/< Tb >• Meanwhile, the MAX/< Tb > 
ratio changes strongly in the range [0.04; 4Q6]K/min and 
noticeably in the range [1; 16]K/min and also depends 
on the K value. Figure 9 shows the results of varying 
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Temperature (K) 


FIG. 6. Comparison between the Tb distribution obtained 
directly from the size distribution used in the simulation and 
the derivative d(ZFC-FC)/dT of the simulated curves. 



FIG. 7. Values of Tb mean, MAX and IP of the simulated 
curves as a function of the scale parameter a for 4.5nm mean 
radius, 16kJ/m^ anisotropy constant and 4:Klmin heating 
rate. 


the heating rate for Rm = lOnm, K = \6kJjw? and 
Ms — T^lkAjm. 

Figure 10 shows a parabolic fit over the IPj < Tb > 
values obtained for all the simulations. The curve is 
universal with small fluctuations due to numeric reso¬ 
lution. The obtained polynomial with fitting errors is 
^4^(cr) = 1.00(2) - 0.21(2)cr - 0.79(2)0-2. 


B. Experimental application 

The Micha’s analysis was conducted on ZFC-FC mea¬ 
surements of a FF of magnetite MNPs suspended in hex¬ 
ane with a concentration of 12{l)g/L. TEM images were 



Scale parameter a 


FIG. 8. Values of MAX (circles) and IP (squares) of the 
simulated curves divided by Tb mean for different MNP radii. 
The behavior is the same for all sizes. 



FIG. 9. MAX (circles) and IP (squares) relative to < Tb > 
for lOnm NPM mean radius at different temperature rates. 


taken in order to determine the size distribution of the 
particles (fig 11). A narrow log-normal number diame¬ 
ter distribution {LnN{x)) was obtained with a 9.54nm 
mean and a 1.73nm standard deviation. The relative 
TEM volume distribution was obtained from this results 
and fitted with a xLnN{x) function obtaining a scale 
parameter a = 0.55(2). 

The ZEC-EC routine was carried with a 2AKjmin rate 
and a SkAjm field on an encapsulated EE sample frozen 
under a IT field in order to obtain an ordered system 
with all the MNP easy axes oriented parallel to the field. 
The ZEC-EC derivative was calculated and fitted with a 
xLnN{x) distribution using the TEM a as a fixed pa¬ 
rameter with a very good correspondence (figure 12). 
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FIG. 10. Parabolic fitting of the universal curve IP/< Tb >• 
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FIG. 11. Size distribution from TEM images. Inset: TEM 
image example with a magnification showing the crystallinity 
of the particles. 


EIG. 12. Log-normal fit on the d(ZEG-EG)/dT derivative. 
Inset: ZEG and EG experimental curves. 


value as the one that maximizes the coincidence between 
TEM and ZFC-FC distributions. In this case, a value of 
34(2)was obtained with a very good correspon¬ 
dence between TEM and ZFC-FC data. This calculation 
implies some approximations: Mg is considered indepen¬ 
dent from the temperature in the region of interest, and 
the relaxation time expression used for the blocking con¬ 
dition 12 considerate only the inversions in the direction 
of the field. While the first approximation is very reason¬ 
ably, the blocking condition expression is accurate only 
in experiments with high jiH/ikT) ratios, where the re¬ 
versal frequency are much smaller for the inversions to 
the antiparallel state. 

Additionally, the IP/< Tb > ratio was calculated ob¬ 
taining a value of 0.7(1), compatible with polynomial ex¬ 
pression obtained from the simulations. 


Figure 13 shows the comparison between Tb distribu¬ 
tion obtained from TEM information and the ZFC-FC 
derivative curve. The translation from TEM volume to 
Tb was made considering the blocking condition in which 
the inversion time of the MNPs is approximately equal 
to the measurement time of the magnetization value: 

t{K, V, h, T) = To exp 

KV{l-hf 
^ fclog(rm/ro)’ 

where tq = l//o is the inverse of the intrinsic inver¬ 
sion frequency. For a known volume distribution, this 
comparison can be used to determine the effective K 


IV. DISCUSSION AND CONCLUSIONS 

The validity of the Micha’s method to determine the 
Tb distribution of non interacting MNPs assembly was 
demonstrated by numerical simulations and experimental 
data analysis. 

A Stoner-Wolfarth model with thermal agitation was 
developed in order to simulate the ZFC-FC curves of 
polysized MNPs assembles. From this simulation it was 
clearly demonstrated that the temperature derivative of 
the ZFC-FC difference is in full coincidence with the Tb 
distribution of the sample, calculated as the inflection 
points of each size ZFC curve. Additionally, it came 
clear from the results that the maximum (MAX) and 
the inflection point (IP) of the polysized ZFC curve are 
affected not only by the mean size of the particles, but by 
the size dispersion. Thus neither IP or MAX are direct 




























































































7 



0 20 40 60 80 


Temperature (K) 

FIG. 13. d(ZFC-FC)/dT derivative together with Tb distri¬ 
bution from TEM volume obtained by fitting K for maximum 
coincidence. 


estimators of the mean Tb- This is an interesting result 
since these values are commonly used in magnetic charac¬ 
terizations and can lead to estimate Tb values far from 
the mean. As an example, for a sample with a = 0.5 
and a heating rate of ^Kjmin^ IP= 0.7 < Tb > and 
MAX= 2.12 < Tb >. Nevertheless, it was found that 
the IP/ < Tb > ratio depends exclusively on cr, while 
MAX/ < Tb > changes with K and the heating-cooling 
rate. This behavior reveals a connection between IP, one 
of the most commonly used Tb criteria and the actual 
mean value of the blocking temperature. For a sample 
with known a, the mean blocking temperature could be 
obtained from the universal curve presented in this work. 
This approach of obtaining size distribution information 
from an universal curve was presented before by Hansen 
et using a rougher model. 

In the present development status, the ZFC-FC simu¬ 
lation algorithm does not include a “measurement time” 
parameter. Just the heating rate is used, so the authors 
assume that the simulated data points represent the val¬ 
ues for “instantaneous” measurements. Therefore, the 
magnetization values obtained from the simulation won’t 
be equivalent to the ones obtained in an experiment with 
the same parameters. 

The Micha’s method was then applied to character¬ 
ize a sample of magnetite nanoparticles coated with oleic 
acid and suspended in hexane. The volume distribution 
of the sample was obtained from TEM analysis show¬ 
ing a narrow log-normal shape with a mean diameter of 
9.54nm and a standard deviation of 1.73nm. In order to 
obtain an ordered system, the ferrofluid was frozen under 
a IT magnetic field. Then, a ZFC-FC routine was car¬ 
ried and a Tg distribution was obtained from the data. 
This distribution was fitted with a xLnN[x\ function us¬ 
ing the scale parameter a obtained from the TEM data 
as a fixed fitting parameter. The high goodness of the 


fitting supports the validity of Micha’s method and the 
low influence of magnetic interaction between particles 
which is consistent with the particle to particle distance 
imposed by the FF concentration. Additionally, the re¬ 
sultant IP/< Tb > values is consistent with the universal 
curve obtained from the simulations. 

Finally, the effective anisotropy constant of the parti¬ 
cles was estimated as the value which gives the maximum 
coincidence between the ZFC-FC Tb distribution and the 
one obtained from the TEM volume. 

The results obtained in this work constitute just a first 
example of the potential of the presented model in combi¬ 
nation with experimental characterization. There is work 
in progress to enhance the simulation algorithm in order 
to include the measurement time as a parameter and to 
considerate both inversions processes in the blocking con¬ 
dition. 

^Giinter Reiss and Andreas Hiitten. Magnetic nanoparticles: Ap¬ 
plications beyond data storage. Nature Materials , 4(10):725-726, 
oct 2005. 

^Quentin A Pankhurst, J Connolly, SK Jones, and JJ Dobson. 
Applications of magnetic nanoparticles in biomedicine. Journal 
of physics D: Applied physics , 36(13):R167, 2003. 

^QA Pankhurst, NTK Thanh, SK Jones, and J Dobson. Progress 
in applications of magnetic nanoparticles in biomedicine. Journal 
of Physics D: Applied Physics , 42(22):224001, 2009. 

"^Ronald E Rosensweig. Heating magnetic fluid with alternating 
magnetic field. Journal of magnetism and magnetic materials, 
252:370-374, 200^ 

^Magdalena Woihska, Jacek Szczytko, Andrzej Majhofer, Jacek 
Gosk, Konrad Dziatkowski, and Andrzej Twardowski. Magnetic 
interactions in an ensemble of cubic nanoparticles: A monte carlo 
study. Physical Review B , 88(14):144421, 2013. 

®Edmund C Stoner and EP Wohlfarth. A mechanism of magnetic 
hysteresis in heterogeneous alloys. Philosophical Transactions 
of the Royal Society of London A: Mathematical, Physical and 
Engineering Sciences , 240(826):599-642, 1948. 

^CP Bean and JD Livingston. Superparamagnetism. Journal of 
Applied Physics , 30(4):S120-S129, 1959. 

^Carolin Schmitz-Antoniak. X-ray absorption spectroscopy on 
magnetic nanoscale systems for modern applications. Reports 
on Progress in Physics , 78(6):062501, 2015. 

® Patricia de la Presa, Yurena Luengo, Victor Velasco, 
MP Morales, M Iglesias, Sabino Veintemillas-Verdaguer, Patri¬ 
cia Crespo, and Antonio Hernando. Particle interactions in liq¬ 
uid magnetic colloids by zero field cooled measurements: Ef¬ 
fects on heating efficiency. The Journal of Physical Chemistry 
C, 119(20):11022-11030, 20lT 

Sandra Sankar, AE Berkowitz, D Dender, JA Borchers, RW Er¬ 
win, SR Kline, and David J Smith. Magnetic correlations in 
non-percolated co-sio 2 granular films. Journal of magnetism 
and magnetic materials , 221(l):l-9, 2000. 

^^F Tournus and A Tamion. Magnetic susceptibility curves of a 
nanoparticle assembly ii. simulation and analysis of zfc/fc curves 
in the case of a magnetic anisotropy energy distribution. Journal 
of Magnetism and Magnetic Materials , 323(9):1118-1127, 2011. 
^^JS Micha, B Dieny, JR Regnard, JF Jacquot, and J Sort. Es¬ 
timation of the CO nanoparticles size by magnetic measurements 
in co/sio 2 discontinuous multilayers. Journal of Magnetism and 
Magnetic Materials , 272:E967-E968, 2004. 

^^H Mamiya, M Ohnuma, I Nakatani, and T Furubayashim. Ex¬ 
traction of blocking temperature distribution from zero-field- 
cooled and field-cooled magnetization curves. IEEE transactions 
on magnetics , 41(10):3394-3396, 2005. 

^"^Jing Ju Lu, Huei Li Huang, and Ivo Klik. Field orientations 





























































and sweep rate effects on magnetic switching of stoner-wohlfarth 
particles. Journal of applied physics , 76(3):1726-1732, 1994. 

Usov and Yu B Grebenshchikov. Hysteresis loops of 
an assembly of superparamagnetic nanoparticles with uniaxial 
anisotropy. Journal of Applied Physics , 106(2):023917, 2009. 
Julian Carrey, Boubker Mehdaoui, and Marc Respaud. Sim¬ 
ple models for dynamic hysteresis loop calculations of magnetic 
single-domain nanoparticles: Application to magnetic hyperther¬ 
mia optimization. Journal of Applied Physics , 109(8):083921, 
2011. 

^^NA Usov. Numerical simulation of field-cooled and zero field- 
cooled processes for assembly of superparamagnetic nanopar¬ 
ticles with uniaxial anisotropy. Journal of Applied Physics , 
109(2):023913, 2011. 

Roger Balian and Dirk Haar. From microphysics to macrophysics: 


methods and applications of statistical physics , volume 2. 
Springer Science &: Business Media, 2006. 

Neel. Influence des fluctuations thermiques a 1 aimantation des 
particules ferromagnetiques. Comptes Rendus de TAcademie des 
Sciences , 228:664-668, 1949. 

^^Amikam Aharoni. Introduction to the Theory of 

Ferromagnetism, volume 109. Oxford University Press, 
2000. 

Lawrence F Shampine and Mark W Reichelt. The matlab ode 
suite. SIAM journal on scientific computing , 18(l):l-22, 1997. 
^^Mikkel Fougt Hansen and Steen Mprup. Estimation of block¬ 
ing temperatures from zfc/fc curves. Journal of Magnetism and 
Magnetic Materials , 203(1):214-216, 1999. 



